Surface soil moisture (SSM) is a crucial parameter for understanding the hydrological process of our earth surface. Passive microwave (PM) technique has long been the primary choice for estimating SSM at satellite remote sensing scales, while on the other hand, the coarse resolution (usually >~10 km) of PM observations hampers its applications at finer scales. Although quantitative studies have been proposed for downscaling satellite PM-based SSM, very few products have been available to public that meet the qualification of 1-km resolution and daily revisit cycles under all-weather conditions. In this study, therefore, we have developed one such SSM product in China with all these characteristics. The product was generated through downscaling of AMSR-E and AMSR-2 based SSM at 36-km, covering all on-orbit time of the two radiometers during 2003-2019. MODIS optical reflectance data and daily thermal infrared land surface temperature (LST) that have been gap-filled for cloudy conditions were the primary data inputs of the downscaling model, in order to achieve the “all-weather” quality for the SSM downscaling outcome. Daily images from this developed SSM product have achieved quasi-complete coverage over the country during April-September. For other months, the national coverage percentage of the developed product is also greatly improved against the original daily PM observations. We evaluated the product against in situ soil moisture measurements from over 2000 professional meteorological and soil moisture observation stations, and found the accuracy of the product is stable for all weathers from clear sky to cloudy conditions, with station averages of the unbiased RMSE ranging from 0.053 vol to 0.056 vol. Moreover, the evaluation results also show that the developed product distinctly outperforms the widely known SMAP-Sentinel (Active-Passive microwave) combined SSM product at 1-km resolution. This indicates potential important benefits that can be brought by our developed product, on improvement of futural investigations related to hydrological processes, agricultural industry, water resource and environment management.
SONG Peilin, ZHANG Yongqiang
Surface melting is the primary reason that affects the mass balance of Greenland ice sheet. At the same time, ice and snow have high albedo, and ice sheet surface melting will cause the difference of radiation energy budget, and then affects the energy exchange between sea-land-air. The high-resolution ice sheet surface melting product provides important information support for the study of Greenland ice sheet surface melting and its response to global climate change. This dataset combined microwave radiometer product and optical albedo product, the daily, winter (June-August) averages and July averages of the former are used for layer-stacking, then Gram-Schmidt Spectral Sharpening was adapted to fuse the layer-stacking results with MODIS GLASS albedo product. The spatial resolution of fusion-results has been downscaled from 25 km to 0.05˚. By employing a threshold-based melt detection approach for each fusion-results pixel, Greenland ice sheet surface melt daily product for 1985, 2000, 2015 (DSSMIS) was generated. The spatial resolution of DSSMIS is higher than that of published data sets at home and abroad. Combined with the advantages of radiometer and albedo data, the spatial details characteristics are enhanced and consistent with the extraction range of the original radiometer products, effectively reducing the noise of the radiometer. DSSMIS’s data type is integer, where 1 is melted, 0 is not melted, 255 is masked area besides Greenland ice sheet, and the data set is stored as *.nc.
WEI Siyi,
This dataset consists of four files including (1) Lake ice thickness of 16 large lakes measured by satellite altimeters for 1992-2019 (Altimetric LIT for 16 large lakes.xlsx); (2) Daily lake ice thickness and lake surface snow depth of 1,313 lakes with an area > 50 km2 in the Northern Hemisphere modeled by a one-dimensional remote sensing lake ice model for 2003-2018 (in NetCDF format); (3) Future lake ice thickness and surface snow depth for 2071-2099 modeled by the lake ice model with a modified ice growth module (table S1.xlsx); (4) A lookup table containing lake IDs, names, locations, and areas. This daily lake ice and snow thickness dataset could provide a benchmark for the estimation of global lake ice and snow mass, thereby improving our understanding of the ecological and economical significance of freshwater ice as well as its response to climate change.
LI Xingdong, LONG Di, HUANG Qi, ZHAO Fanyu
The freeze/thaw status of the near-surface soil is the water-ice phase transition that occurred at the top soil layer. It is an important indicator as a giant on-off “switch” of the land surface processes including water, energy, and carbon exchanges between the land surface and atmosphere. The freeze/thaw status is an essential variable for understanding how the ecosystem responds to and affects global changes. This dataset is based on the AMSR-E and AMSR2 passive microwave brightness temperature data, and the freeze-thaw discriminant function algorithm is used to generate the global near-surface soil freeze-thaw status with a spatial resolution of grids at 0.25° from 2002 to 2019. The dataset can be used for the analysis of the spatial distribution and trend changes of global freeze-thaw cycles, such as the freeze/thaw onset dates and duration. It provides data support for understanding the interaction mechanism between the land surface freeze-thaw cycle and the land-atmosphere exchanges under the context of global changes.
Zhao Tianjie
The freeze/thaw status of the near-surface soil is the water-ice phase transition that occurred at the top soil layer. It is an important indicator as a giant on-off “switch” of the land surface processes including water, energy, and carbon exchanges between the land surface and atmosphere. The freeze/thaw status is an essential variable for understanding how the ecosystem responds to and affects global changes. This dataset is based on the AMSR-E, AMSR2 passive microwave brightness temperature data and MODIS optical remote sensing data. The freeze-thaw discriminant function algorithm and downscaling algorithm are used to generate the global near-surface soil freeze-thaw status with a spatial resolution of grids at 0.05° from 2002 to 2017. The dataset can be used for the analysis of the spatial distribution and trend changes of global freeze-thaw cycles, such as the freeze/thaw onset dates and duration. It provides data support for understanding the interaction mechanism between the land surface freeze-thaw cycle and the land-atmosphere exchanges under the context of global changes.
Zhao Tianjie, ZHANG Ziqian
There are many lakes in the Qinghai Tibet Plateau. The glacial phenology and duration of lakes in this region are very sensitive to regional and global climate change, so they are used as the key indicators of climate change research, especially the comparative study of the three polar environmental changes of the earth. However, due to its poor natural environment and sparse population, there is a lack of conventional field measurement of lake ice phenology. The lake ice was monitored with a resolution of 500 meters by using the normalized difference snow index (NDSI) data of MODIS. The traditional snow map algorithm is used to detect the lake daily ice amount and coverage under the condition of sunny days, and the lake daily ice amount and coverage under the condition of cloud cover are re determined through a series of steps based on the spatiotemporal continuity of the lake surface conditions. Through time series analysis, 308 lakes larger than 3km2 are identified as effective records of lake ice range and coverage, forming a daily lake ice range and coverage data set, including 216 lakes.
QIU Yubao
The Antarctic ice sheet elevation data were generated from radar altimeter data (Envisat RA-2) and lidar data (ICESat/GLAS). To improve the accuracy of the ICESat/GLAS data, five different quality control indicators were used to process the GLAS data, filtering out 8.36% unqualified data. These five quality control indicators were used to eliminate satellite location error, atmospheric forward scattering, saturation and cloud effects. At the same time, dry and wet tropospheric, correction, solid tide and extreme tide corrections were performed on the Envisat RA-2 data. For the two different elevation data, an elevation relative correction method based on the geometric intersection of Envisat RA-2 and GLAS data spot footprints was proposed, which was used to analyze the point pairs of GLAS footprints and Envisat RA-2 data center points, establish the correlation between the height difference of these intersection points (GLAS-RA-2) and the roughness of the terrain relief, and perform the relative correction of the Envisat RA-2 data to the point pairs with stable correlation. By analyzing the altimetry density in different areas of the Antarctic ice sheet, the final DEM resolution was determined to be 1000 meters. Considering the differences between the Prydz Bay and the inland regions of the Antarctic, the Antarctic ice sheet was divided into 16 sections. The best interpolation model and parameters were determined by semivariogram analysis, and the Antarctic ice sheet elevation data with a resolution of 1000 meters were generated by the Kriging interpolation method. The new Antarctic DEM was verified by two kinds of airborne lidar data and GPS data measured by multiple Antarctic expeditions of China. The results showed that the differences between the new DEM and the measured data ranged from 3.21 to 27.84 meters, and the error distribution was closely related to the slope.
HUANG Huabin
This biophysical permafrost zonation map was produced using a rule-based GIS model that integrated a new permafrost extent, climate conditions, vegetation structure, soil and topographic conditions, as well as a yedoma map. Different from the previous maps, permafrost in this map is classified into five types: climate-driven, climate-driven/ecosystem-modified, climate-driven/ecosystem protected, ecosystem-driven, and ecosystem-protected. Excluding glaciers and lakes, the areas of these five types in the Northern Hemisphere are 3.66×106 km2, 8.06×106 km2, 0.62×106 km2, 5.79×106 km2, and 1.63×106 km2, respectively. 81% of the permafrost regions in the Northern Hemisphere are modified, driven, or protected by ecosystems, indicating the dominant role of ecosystems in permafrost stability in the Northern Hemisphere. Permafrost driven solely by climate occupies 19% of permafrost regions, mainly in High Arctic and high mountains areas, such as the Qinghai-Tibet Plateau.
RAN Youhua, M. Torre Jorgenson, LI Xin, JIN Huijun, Wu Tonghua, Li Ren, CHENG Guodong
The “Long-term series of daily global snow depth” was produced using the passive microwave remote sensing data. The temporal range is 1980~2018, and the coverage is the global land. The spatial resolutions is 25,067.53 m and the temporal resolution is daily. A dynamic brightness temperature gradient algorithm was used to derive snow depth. In this algorithm, the spatial and temporal variations of snow characteristics were considered and the spatial and seasonal dynamic relationships between the temperature difference between 18 GHz and 36 GHz and the measured snow depth were established. The long-term sequence of satellite-borne passive microwave brightness temperature data used to derive snow depth came from three sensors (SMMR, SSM/I and SSMI/S), and there is a certain system inconsistency among them. So, the inter-sensor calibration was performed to improve the temporal consistency of these brightness temperature data before snow depth derivation. The accuracy analysis shows that the relative deviation of Eurasia snow depth data is within 30%. The data are stored as a txt file every day, each file is a 1383*586 snow depth matrix, and each snow depth represents a 25,067.53m* 25,067.53m grid. The projection of this data is EASE-Grid, and following is the file header which describes the projection detail. File header: ncols 1383 nrows 586 xllcorner -17334193.54 yllcorner -7344787.75 cellsize 25,067.53 NODATA_value -1
CHE Tao, LI Xin, DAI Liyun
This dataset was derived from long-term daily snow depth in China based on the boundary of the three-river-source area. The snow depth ranges from 0 to 100 cm, and the temporal coverage is from January 1 1980 to December 31 2020. The spatial and temporal resolutions are 0.25o and daily, respectively. Snow depth was produced from satellite passive microwave remote sensing data which came from three different sensors that are SMMR, SSM/I and SSMI/S. Considering the systematic bias among these sensors, the inter-sensor calibrations were performed to obtain temporal consistent passive microwave remote sensing data. And the long-term daily snow depth in China were produced from this consistent data based on the spectral gradient method.For header file information, refer to the data set header.txt.
DAI Liyun
This data set uses SMMR (1979-1987), SSM / I (1987-2009) and ssmis (2009-2015) daily brightness temperature data, which is generated by double index (TB V, SG) freeze-thaw discrimination algorithm. The classification results include four types: frozen surface, melted surface, desert and water body. The data covers the source area of three rivers, with a spatial resolution of 25.067525 km. It is stored in geotif format in the form of ease grid projection. Pixel values represent the state of freezing and thawing: 1 for freezing, 2 for thawing, 3 for deserts, 4 for water bodies. Because all TIF files in the dataset describe the scope of Sanjiangyuan National Park, the row and column number information of these files is unchanged, and the excerpt is as follows (where the unit of cellsize is m): ncols 52 nrows 28 cellsize 25067.525 nodata_value 0
River lake ice phenology is sensitive to climate change and is an important indicator of climate change. 308 excel file names correspond to Lake numbers. Each excel file contains six columns, including daily ice coverage information of corresponding lakes from July 2002 to June 2018. The attributes of each column are: date, lake water coverage, lake water ice coverage, cloud coverage, lake water coverage and lake ice coverage after cloud treatment. Generally, the ice cover area ratio of 0.1 and 0.9 is used as the basis to distinguish the lake ice phenology. The excel file contained in the data set can further obtain four lake ice phenological parameters: Fus, fue, bus, bue, and 92 lakes. Two parameters, Fus and bue, can be obtained.
QIU Yubao
The global Cryosat-2 GDR dataset is generated by the European Space Agency (ESA); it has a temporal coverage from 2010 to 2016 and covers the globe. On April 8, 2010, the ESA launched the Cryosat-2 high-tilt polar orbit satellite. The satellite is equipped with an SAR Interferometer Radar Altimeter (SIRAL), which is mainly used to monitor polar ice thickness and sea ice thickness changes, and, furthermore, to study the effects of melting polar ice on global sea level rise and that of global climate change on Antarctic ice thickness. The altimeter operates in the Ku-band and at a frequency of 13.575 GHz, it includes three measurement modes. One is a low-resolution altimeter measurement mode (LRM) that points to the subsatellite point to obtain all surface observations for land, sea, and ice sheets; its processing is similar to ENVISAT/RA-2, with an orbital resolution of 5 to 7 km. The second is the Synthetic Aperture Radar (SAR) measurement mode, which is mainly used to improve the accuracy and resolution of sea ice observations; it can make the resolution along the orbit reach approximately 250 m. The third is the Interferometric Synthetic Aperture Radar (InSAR), which is mainly used to improve the accuracy of areas with complex terrain such as the edges of ice sheets or ice shelves. The CryoSat -2/SIRAL data products mainly include 0-level data, 1b-level data, 2-level data and high-level data. The Cryosat-2/SIRAL products consist of two files: an XML head file (.HDR) and a data product file (.DBL). The HDR file is an auxiliary ASCII file for fast identification and retrieval of the data files. 1b-level products are stored separately according to the measurement modes, and the data recording formats of different modes are also different. Each waveform in LRM mode and SAR mode has 128 sampling points, while that in SARIn mode has 512 sampling points. 2-level GDR products are available for most scientific applications, including measurement time, geographic location, altitude, and more. In addition, the altitude information in GDR products has been obtained through instrumental calibration, transmission delay corrections, geometric corrections, and geophysical corrections (such as atmospheric corrections and tidal corrections). The GDR products are single global full-track data, that is, the measurement results of the three modes. After different processing, they are combined in chronological order; thereby, the data recording formats are unified. The data in the three modes use different waveform retracking algorithms to obtain altitude values. In the latest updated Baseline C data, the LRM mode data use three algorithms: Refined CFI, UCL and Refined OCOG.
SHEN Guozhuang, FU Wenxue
This dataset includes passive microwave remote sensing brightness temperatures data for longitude and latitude projections and 0.25 degree resolution from 2002 to 2008 in China. 1. Data processing process: NSIDC produces AMSR-E gridded brightness temperature data by interpolating AMSR-E data (6.9 GHz, 10.7 GHz, 18.7 GHz, 23.8 GHz, 36.5 GHz, and 89.0 GHz) to the output grids from swath space using an Inverse Distance Squared (ID2) method. 2. Data format: Brightness temperature files: two-byte unsigned integers, little-endian byte order Time files: two-byte signed integers, little-endian byte order 3. Data naming: ID2rx-AMSRE-aayyyydddp.vnn.ccc (China-ID2r1-AMSRE-D.252002170A.v03.06V) ID2 Inverse Distance Squared r1 Resolution 1 swath input data AMSRE Identifies this an AMSR-E file D.25 Identifies this as a quarter degree file yyyy Four-digit year ddd Three-digit day of year p Pass direction (A = ascending, D = descending) vnn Gridded data version number (for example, v01, v02, v03) ccc AMSR-E channel indicator: numeric frequency (06, 10, 18, 23, 36, or 89) followed by polarization (H or V) 4. Cutting range: Corner Coordinates: Upper Left (60.0000000, 55.0000000) (60d 0'0.00 "E, 55d 0'0.00" N) Lower Left (60.0000000, 15.0000000) (60d 0'0.00 "E, 15d 0'0.00" N) Upper Right (140.0000000, 55.0000000) (140d 0'0.00 "E, 55d 0'0.00" N) Lower Right (140.0000000, 15.0000000) (140d 0'0.00 "E, 15d 0'0.00" N) Center (100.0000000, 35.0000000) (100d 0'0.00 "E, 35d 0'0.00" N) Origin = (60.000000000000000, 55.000000000000000) 5. Data projection: GEOGCS ["WGS 84", DATUM ["WGS_1984", SPHEROID ["WGS 84", 6378137,298.257223563, AUTHORITY ["EPSG", "7030"]], TOWGS84 [0,0,0,0,0,0,0], AUTHORITY ["EPSG", "6326"]], PRIMEM ["Greenwich", 0, AUTHORITY ["EPSG", "8901"]], UNIT ["degree", 0.0174532925199433, AUTHORITY ["EPSG", "9108"]], AUTHORITY ["EPSG", "4326"]]
Mary Jo Brodzik, Matthew Savoie, Richard Armstrong, Ken Knowles
Snow is a significant component of the ecosystem and water resources in high-mountain Asia (HMA). Therefore, accurate, continuous, and long-term snow monitoring is indispensable for the water resources management and economic development. The present study improves the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua satellites 8 d (“d” denotes “day”) composite snow cover Collection 6 (C6) products, named MOD10A2.006 (Terra) and MYD10A2.006 (Aqua), for HMA with a multistep approach. The primary purpose of this study was to reduce uncertainty in the Terra–Aqua MODIS snow cover products and generate a combined snow cover product. For reducing underestimation mainly caused by cloud cover, we used seasonal, temporal, and spatial filters. For reducing overestimation caused by MODIS sensors, we combined Terra and Aqua MODIS snow cover products, considering snow only if a pixel represents snow in both the products; otherwise it is classified as no snow, unlike some previous studies which consider snow if any of the Terra or Aqua product identifies snow. Our methodology generates a new product which removes a significant amount of uncertainty in Terra and Aqua MODIS 8 d composite C6 products comprising 46 % overestimation and 3.66 % underestimation, mainly caused by sensor limitations and cloud cover, respectively. The results were validated using Landsat 8 data, both for winter and summer at 20 well-distributed sites in the study area. Our validated adopted methodology improved accuracy by 10 % on average, compared to Landsat data. The final product covers the period from 2002 to 2018, comprising a combination of snow and glaciers created by merging Randolph Glacier Inventory version 6.0 (RGI 6.0) separated as debris-covered and debris-free with the final snow product MOYDGL06*. We have processed approximately 746 images of both Terra and Aqua MODIS snow containing approximately 100 000 satellite individual images. Furthermore, this product can serve as a valuable input dataset for hydrological and glaciological modelling to assess the melt contribution of snow-covered areas. The data, which can be used in various climatological and water-related studies, are available for end users at https://doi.org/10.1594/PANGAEA.901821 (Muhammad and Thapa, 2019).
SHER Muhammad
The aerosol optical thickness data of the Arctic Alaska station is based on the observation data products of the atmospheric radiation observation plan of the U.S. Department of energy at the Arctic Alaska station. The data coverage time is updated from 2017 to 2019, with the time resolution of hour by hour. The coverage site is the northern Alaska station, with the longitude and latitude coordinates of (71 ° 19 ′ 22.8 ″ n, 156 ° 36 ′ 32.4 ″ w). The source of the observed data is retrieved from the radiation data observed by mfrsr instrument. The characteristic variable is aerosol optical thickness, and the error range of the observed inversion is about 15%. The data format is NC format. The aerosol optical thickness data of Qomolangma station and Namuco station in the Qinghai Tibet Plateau is based on the observation data products of Qomolangma station and Namuco station from the atmospheric radiation view of the Institute of Qinghai Tibet Plateau of the Chinese Academy of Sciences. The data coverage time is from 2017 to 2019, the time resolution is hour by hour, the coverage sites are Qomolangma station and Namuco station, the longitude and latitude coordinates are (Qomolangma station: 28.365n, 86.948e, Namuco station Mucuo station: 30.7725n, 90.9626e). The source of the observed data is retrieved from the radiation data observed by mfrsr instrument. The characteristic variable is aerosol optical thickness, and the error range of the observed inversion is about 15%. The data format is TXT.
QIU Yubao
Using the Modis1B data of 11 scenes from 2003 to 2013 (the ice shelf Modis1B data published on the NSIDC website), the surface velocity of the Antarctic Amery Ice Shelf was extracted by the subpixel cross-correlation method, the ice velocity was extracted by the COSI-Corr software, and then the time sequence of annual average velocities for nearly ten years was obtained. Due to the lack of field observations in the study area, the accuracy of the ice flow results was estimated by using the offset value of the stable region, and the ice flow error was approximately ±50 m/year. The ice velocity data date from 2003 to 2013, the temporal resolution is one year, and the data cover the Amery area with a spatial resolution of 500 m. A GeoTIFF file of velocity data is stored every year. For details regarding the data, please refer to the Amery Ice Flow Field - Data Description.
JIANG Liming
As the main parameter in the land surface energy balance, surface temperature indicates the degree of land-atmosphere energy and water transfer and is widely used in research on climatology, hydrology and ecology. In the study of frozen soil, climate is one of the decisive factors for the existence and development of frozen soil. The surface temperature is the main climatic factor affecting the distribution of frozen soil and affects the occurrence, development and distribution of frozen soil. It is the upper boundary condition for modelling frozen soil and is significant to the study of hydrological processes in cold regions. The data set was based on the DEM and observation station data of the Tibetan Plateau Engineering Corridor and analysed the changing trend of surface temperature on the Tibetan Plateau from 2000 to 2014. Using the surface temperature data products MOD11A1/A2 and MYD11A1/A2 of MODIS aboard Terra and Aqua, the surface temperature information under cloud cover was reconstructed based on the spatio-temporal information of the images. The reconstruction information and surface temperature representativeness problems were analysed using information obtained from 8 sites, including the Kunlun Mountains (wetland, grassland), Beiluhe (grassland, meadow), Kaixinling (meadow, grassland), and Tanggula Mountain (meadow, wetland). According to the correlation coefficient (R2), root-mean-square error (RMSE), mean absolute error (MAE) and mean deviation (MBE), the following results were obtained: (1) the reconstruction accuracy of MODIS surface temperature under cloud cover is higher when it is based on spatio-temporal information; (2) the weighted average representation is the best when generalizing four observations of Terra and Aqua. By analysing the reconstruction of MODIS surface temperature information and representativeness problems, the average annual MODIS surface temperature data of the Tibetan Plateau and the engineering corridor from 2000 to 2010 were obtained. According to the data set, the surface temperature from 2000 to 2010 also experienced volatile rising trends from 2000 to 2010, which is basically consistent with the changing trend of the climate change in the permafrost regions of the Tibetan Plateau and the Qinghai-Tibet Engineering Corridor.
NIU Fujun, YIN Guoan
The variation in the duration of snow on the Tibetan Plateau is relatively great, and the high mountainous areas around the plateau are rich in snow and ice resources. Taking full account of the terrain of the Tibetan Plateau and the snow characteristics in the mountains, the data set adopted AVHRR data to gradually realize generating data products for daily, ten-day, and monthly snow cover areas while maintaining the snow classification accuracy. These data included the daily/10-day/monthly snow cover area data for the Tibetan Plateau from 2007 to 2015, the average accuracy of which is 0.92. It can provide reliable data for snow changes during the historical periods of the Tibetan Plateau.
QIU Yubao
The proportion data set of daily cloudless MODIS snow cover area in babaohe river basin (2008.1.1-2014.6.1) was obtained after cloud removal processing using a cloud removal algorithm based on cubic spline function interpolation on the basis of daily cloudless MODIS snow cover product-mod10a1 (tang zhiguang, 2013). This data set adopts the projection method of UTM (horizontal axis isometric cutting cylinder), with a spatial resolution of 500m, and provides Daily Snow Albedo daily-sad results for the babao river basin.The data set is a daily file from January 1, 2008 to June 1, 2014.Each file is the snow albedo result of the day, with a value of 0-100 (%), is the ENVI standard file, and the naming rule is: mod10a1.ayyyyddd_h25v05_snow_sad_grid_2d_reproj_babaohe_nocloud.img, where YYYY represents the year, DDD stands for Julian day (001-365/366).The file can be opened directly with ENVI or ARCMAP software. The original MODIS snow cover data products processed by declouding are derived from MOD10A1 products processed by the us national snow and ice data center (NSIDC). This data set is in HDF format and USES sinusoidal projection. The attributes of the cloud-free MODIS albedo data set (2008.1.1-2014.1.1) in babaohe river basin are composed of the spatial and temporal resolution, projection information and data format of the dataset.
WANG Jian, PAN Haizhu
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS 0931-4967287 poles@itpcas.ac.cnLinks
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved | No.11010502040845
Tech Support: westdc.cn